Electron-nuclear dynamics in a quantum dot under non-unitary electron control 
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We introduce a method for solving tlie problem of an externally controlled electron spin in a 
quantum dot interacting with host nuclei via the hyperfine interaction. Our method accounts 
for generalized (non-unitary) evolution effected by external controls and the environment, such 
as coherent lasers combined with spontaneous emission. As a concrete example, we develop the 
microscopic theory of the dynamics of nuclear-induced frequency focusing as first measured in Science 
317, 1896 (2007); we find that the nuclear relaxation rates are several orders of magnitude faster 
than those quoted in that work. 



The nuclear environment in III-V quantum dots has 
been recognized in recent years as the main source 
of decoherence for the electron spin and thus consti- 
tutes an important hurdle for quantum technologies with 
these systems. The microscopic dynamics of the closed 
electron-nuclear spin system have been investigated in 
important recent theoretical contributions [U [2] . In these 
works [1], controls have been represented as ideal, unitary 
rotations of the electron spin, and the nuclear polariza- 
tion along the external field is taken to be unaltered dur- 
ing the electron evolution. This no longer is the case in 
experiments involving controls that couple the system to 
an additional bath, which can exchange polarization with 
the system. Such experiments are relevant because inco- 
herent interactions are needed to initialize and read out 
the system. These experiments in quantum dots (QDs) 
observed dynamic nuclear polarization and nuclear feed- 
back effects [SHS]. While the details of the various ex- 
periments differ, the main common feature is that an 
external control interacts with the electron, and through 
the hyperfine interaction the nuclear spins are also par- 
tially polarized. The theories employed to describe such 
experiments are usually in the form of rate equations 
and some sort of Fermi golden rule, and typically invoke 
phenomenological terms. Other theories iij employed a 
more microscopic approach, but without including, e.g., 
feedback and a complete treatment of control fields. 

In this Letter, we develop a theory that addresses such 
experiments involving non-unitary evolution of the elec- 
tron while still treating the electron-nuclear interaction 
microscopically. We make use of the operator sum repre- 
sentation of quantum evolution and its simplified form in 
the spin vector (SV) representation and develop a theory 
that is perturbative with respect to the hyperfine cou- 
pling. We develop both Markovian and non-Markovian 
treatments, and by comparison of the two we establish 
the validity regime of the Markovian approximation. 

In order to illustrate the power of our approach, we 
apply it to the experiment of Ref.(3], a proper micro- 
scopic theory of which is lacking to date. This experi- 
ment demonstrated nuclear-induced focusing of the elec- 
tron precession rates in a QD ensemble through the feed- 



back dynamics of the electron and nuclear spins. This 
mechanism is largely driven by non-unitary evolution of 
the electron spin, making it difficult to solve conven- 
tional Master Equations to analyze the dynamics. In- 
stead, a phenomenological treatment was introduced in 
the Supporting Online Material of [3 and further devel- 
oped in [5]. Our microscopic solution does not invoke 
phenomenological quantities and provides a unified de- 
scription of the experiments in [31 |S] . One of our striking 
results is that the nuclear relaxation process is several 
orders of magnitude faster than what is used in [31 [S] . 

The system we consider is a single electron trapped in 
a QD and subject to an external in-plane static magnetic 
field Bz , which splits the spin states along the z direction. 
The electron interacts through the hyperfine contact in- 
teraction with N nuclear spins in the QD {N « 10^). 
There is also an external time dependent field acting on 
the electron, as well as the photon bath that drives spon- 
taneous emission (electron- hole recombination). The to- 
tal Hamiltonian is H = Hq + H^f + Hp + Hrad, where 

Ho = ujj, + eT\T){T\+iJr,J2fi (1) 

Hhf = 5] A,/2(5+/l+5_/;) (2) 

Hp = n{t)\x){T\+h.c. (3) 

Hrad = J2,9ki\z){T\ + \z){T\)ale^'^-' + h.c. (4) 

In Eqs. (l)-(4), Sj{Ij) is the electron (nuclear) spin op- 
erator along the j axis, S± = Sx ±iSy, fl(t) contains the 
pulse information, |z)(|z)) is the spin up (down) state 
along the B-field direction, \x) is the spin down state 
along the optical axis x, \T) is the excited trion state, g 
is the coupling to the radiation bath, and is the bath 
photon creation operator. In the hyperfine Hamiltonian 
Hhf , the first term is referred to as the Overhauser term, 
while the second is called the 'flip-fiop' term. 

The couplings ^l{t) arise from periodic, ultrafast laser 
fields like those used in [3]. Since Ai/oJe ^ 1 for mod- 
erate magnetic fields, we give a perturbative treatment 
in this small parameter. We first focus on the zeroth 
order solution (electron spin periodically driven without 
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FIG. 1: Color online. Left: Experimental setup of [S) |5]. 
Right: Relevant QD states and polarization selection rules. 



nuclear coupling). The primary effect of the pulses on 
the electron spin dynamics is the creation or destruction 
of spin polarization depending on the spin state. This 
arises from the selection rules of the three-level system 
in conjunction with the perpendicular external magnetic 
field. For concreteness, we consider cr^ pulses, in which 
case only \x) is coupled to the light and excited by it to 
the trion, see Fig. [T] Depending on the pulse parameters, 
a certain population is moved to the trion. This popu- 
lation subsequently decays back to the spin subspace via 
spontaneous emission of a photon. Due to the B-field, 
the population decays equally to the \x) and \x) states, 
changing the electron spin polarization [7]. 

This physics describes non unitary evolution of the 
electron spin due to the coupling of the system to the 
photon bath. To describe this mathematically in the spin 
subspace we need a generalization of the usual unitary 
evolution operator to a set of so-called Kraus operators 
{Ej} which transform the density matrix as EjpEj 
[8j. These can be found by solving for the non-unitary 
part of the evolution of an arbitrary initial system density 
matrix and relating it to the final density matrix. Fol- 
lowing this standard procedure [S] we find the following 
Kraus operators in the \x), \x) basis 
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q = QoC c..^ - V 2 

the probability to go from \x) to \x 



The parameter is 
); it is related to the 
pulse area and takes values from to 1. The quantity 
(1 — ^o) is related to the probability of population re- 
maining in the trion state after the passage of the pulse, 
and thus qo quantifies the deviation from unitary dynam- 
ics in the qubit subspace (for unitary evolution qo ~ 1). 
The parameter (j) is the spin rotation angle caused by the 
pulse and is a function of the detuning. We have there- 
fore found 2-d matrices to describe the more complicated 
dynamics of the pulse followed by spontaneous emission. 

In between pulses the evolution is simply Larmor pre- 
cession under B^, given hy U — ■^^ , where Tr 
is the period of the pulse train. We are interested in 



finding the steady state electron spin. For this, the SV 
representation {Sej = 2Tr{pSj)) is most convenient as 
all the operations act on the left side of the SV. As a 
result of the non-unitarity of the evolution, in addition 
to the transformation of the SV a new contribution is 
generated at each cycle: 



(6) 



where we found that {Ye)ij 



^ 2EfTr \s,EiUSjWE\ 

In the limit n — > oo the 

steady state is S'f"'^ = {l — Y(.)^^Ke (explicit expression 
is in PI). We therefore see that a 3 x 3 matrix. Ye, and 
a three-dimensional vector, K^, are the quantities that 
determine the dynamics of the electron spin. Because its 
structure is convenient we use the equivalent and more 
compact 4x4 matrix that contains all the information: 
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In this 4-d representation, the steady-state SV Se°°^ = 
(1, si'^\ si'^\ si^^) is the eigenvector of 1 — J^e with 
eigenvalue 0. This more compact representation will 
prove very useful when we introduce the nuclear spin. 

Having solved the zeroth order problem, we proceed 
to the inclusion of the hyperfine interaction. For simplic- 
ity we assume that the nuclear spin has / = 1/2. The 
nuclear spins affect each other through their interactions 
with the electron. When Ay/N/uJe ^ 1, where ^ is a 
typical value of Ai, fiip-flops occur slowly so that multi- 
nuclear effects such as dark state saturation [B] are neg- 
ligible, and the primary effect of the nuclear spins on the 
electron is a shift of the precession frequency through the 
Overhauser term (Overhauser shift). Therefore, we con- 
sider first a single nuclear spin interacting with the elec- 
tron and incorporate multinuclear effects by shifting the 
electron Zeeman frequency by an amount proportional to 
the net nuclear polarization ^3, 5J. 

For a single nuclear spin interacting with the electron 
spin via the hyperfine Hamiltonian we use a SV represen- 
tation, which in this case is 15-d. For the type of control 
used in Ref. [3J [S] , there are no nuclear effects during 
the ultrashort (i.e., broadband) pulses, which do not dis- 
tinguish between the electron spin eigenstates along the 
field Bz- Therefore, the Kraus operators are simply ten- 
sor products between the Ej 's of Eq. ^ and the identity. 
Following the same prescription as for the single spin, we 
define a 16-d SV 5, = 4Tr (pd) , where p is the 4 x 4 
density matrix of the two spins and the generators Gi are 
tensor products of spin operators (including the identity) 
Gik+e = Sk ® It, where run from to 3. With our 
conventions, Sq = 1. The 16-d analogue of 3^e is given by 
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In general S is not 
but 



simply a tensor product of the two individual SVs, 
contains quantum correlations (entanglement). 

The pulses are expected to 'interrupt' the electron- 
nuclear evolution only for go <C 1, while entanglement 
will build up when Qo ^ 1. Therefore a Markovian ap- 
proximation should be sufhcient for short pulse train pe- 
riods and pulses of strength ^ (as in [31 [S]) . 

Markovian approximation~To find an effective relax- 
ation rate for the nuclear spin, we use the equation 
S(t + Tji) = yS{t) for the 16-d case. In the Markovian 
approximation, we only keep the separable (tensor prod- 
uct) part of S, i.e., S = (8) where we have used 
that the timescales of evolution for the electron and the 
nuclei are quite different (TU] , so that we can assume that 
the electron steady state is reached quickly compared to 
the nuclear dynamics [31 [51 [TTl [T^]. The equation for 
the 4-d nuclear SV is then Sn{t + Tr) = 3^n5„(t), where 
yn explicitly contains electron SV components. Since 
the nuclear evolution is much slower than the pulse rep- 
etition rate, we can coarse grain this equation, and ob- 
tain a differential equation for the nuclear SV, = 
T^{yn - 1)5„, which gives Sn{t) = e(^"-i)*/^«5„(0). 
For small flip-flop coupling (but keeping the Overhauser 
term to all orders), we find the two smallest eigenvalues 
of 1 — Xi to be Ai = and 



l + gg^ + (5g-l)cos(^) 

a.,M + 5L + (^e^-l)cOs(^) 



, (8) 



where is the length of the electron steady state SV, and 
for brevity we have suppressed the superscript c» [13] . 
The zero eigenvalue corresponds to the nuclear steady- 
state SV, which to leading order in the flip-flop term is 
5, 



(oo) _ 



{1,0,0, Sk°:^^) where 



c(oo) 



1 + 52 z 



l)cos(^) 



(9) 



The nonzero eigenvalue A2 gives the nuclear relaxation 
rate 7„ = A2/T'fl. The single nucleus spin flip rates, 
which are generally different in the presence of nonzero 
polarization [S], are w]. = 7„(1 ± S'i'^'')/2, where w\, 
(w^ ) is the rate to flip from down (up) to up (down) . Fig. 
[2]shows that our rates are orders of magnitude larger than 
those of [31 [S] . The heuristic expressions only took into 
account that the relaxation rates should vanish when 
is a multiple of the electron spin precession period as well 
as the overall scale factor A^/wf. The first of these fea- 
tures arises because an electron spin synchronized with 
the pulses is unaffected by them so that no nuclear re- 
laxation takes place. The scale factor is fixed by noting 
that energy conservation leads to a suppression of hyper- 
fine flip-flops when oj^ ^ ujn', only virtual flip-flops are 
allowed and since these must come in pairs, their effect 
is second-order in A/ajg. 
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FIG. 2: Color online. Plot of the log of the nuclear relaxation 
rates multiplied by Tr from our current theory (upper set of 
curves) with our Markovian approximation (solid) and exact 
numerical (dashed) as function of qo- The lower curves are 
based on Refs. [31 H] and are several orders of magnitude less. 
Pulse parameters are (f) = n/2, Tr = 3900.3/tJe, ^/<^e = 10~^. 



Our theory reveals an additional dependence of the re- 
laxation rate 7„ on the orientation of the electron spin 
that was overlooked by |31|3] . When the electron SV is 
transverse to the B-field {Se^z ~ and 0), flip-flops 

are not suppressed by energy conservation and angular 
momentum is freely transferred from the electron to the 
nuclei, leading to a strong enhancement of 7„. This is 
also clear from Eq. ([S]) where the denominator is close 
to zero when Se,z ~ while the numerator remains fi- 
nite due to Se ^ 0. These conditions are realized in the 
regime most relevant for the experiments in Refs. [31 [5] 
where qo 1, that is, when the pulses drive most of the 
population out of the qubit subspace, re-orienting the 
electron spin along the optical (x) axis. Note that the 
enhancement of 7„ depends crucially on the openness 
of the system since the photon bath acts as an angu- 
lar momentum reservoir. Our predicted timescale can be 
checked experimentally by measuring in a single QD the 
frequency of the pump-probe signal at various timescales. 
By systematically varying the B-field and pulse parame- 
ters the relaxation rates could be mapped out as a func- 
tion of the parameters. 

The probability distribution for the net multinuclear 
polarization m/2 is obtained from a kinetic equation for 
m, which is the difference in the number of spins pointing 
up and down: 



dP{m) 
dt 



E 



w±{m) 



^P(m±2)wzp(m±2) 



P{m) (10) 
'N±m 



1 



where w± (m) are the rates in the presence of nuclear 
polarization to/2. These are found by implementing the 
Overhauser shift, w±{m) = w^.^uJe + mA), where we have 
assumed equal couplings for all nuclear spins [14j. Ex- 
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FIG. 3: Color online. The nuclear polarization probability 
distribution with A'^ = 5000 nuclei and pulse parameters = 
7r/4,(7o = 0.2, A/uje = 10"^ and for (a) Tr = 797.9/cJe and 
(b) Tr = 800.3/aje, from our Markovian theory (blue/solid 
line), the theory from Ref. 3 (red/dotted) and the one from 
Ref . [5] (green / dashed) . 
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FIG. 4: Color online. Nuclear steady state SV z com- 
ponent as function of for Markovian approximation 
Eq. Q (red/dotted), analytical non-Markovian, Eq. (Ill 
(green/solid) and exact numerical (black/dashed) for = 
■k/2,Tr = 3900.3/a;e, A/we = 10"^ Inset shows S„,:,, S„,y. 



amples of the resulting distribution are shown in Fig. [3] 
for typical values of the parameters. In general, large 
Tr results in more peaks in P{m) and thus gives rise 
to a greater degree of nuclear state "narrowing" (T2* is 
enhanced). Furthermore, the sharpest peaks occur at 
values of m such that (wg + mA)Tji is an odd integer 
multiple of tt, and the locations of these peaks can be 
controlled by adjusting tOeTn. A systematic exploration 
of the parameter space can help tailor the nuclear state. 

Beyond the Markovian approximation-Our analysis 
above provides analytic expressions for the nuclear dy- 
namics in the Markovian approximation, an approach 
valid for 1 (see Fig. [2]). Our formalism how- 

ever is not inherently Markovian, and we now present 
an analytical non-Markovian expression for the nuclear 
steady state. We return to the 16-d matrix y and per- 
form a perturbative expansion in the coupling which is a 
controlled approximation in the hyperfine coupling. The 
steady state nuclear SV turns out to be 

Sir^ = c(l + ql- 2qo COS (t>, (1 - 9o) tan(c^eTfl,/2), 
qoiVeTn sin (j)), (11) 



where c is in [TS]. Nonzero x,y components arise from 
expanding the Overhauser interaction in addition to the 
flip-flop in deriving Eqn. (111. Fig. |4] shows that the 



dynamics become less Markovian as the pulses become 
more unitary (go — >■ !)• 

In conclusion, we have developed a formalism for an- 
alyzing experiments with generalized, non-unitary con- 
trols on the electron spin confined in a QD and coupled 
to the host nuclei. By applying it to the experiments of 
[31 |S] we have found that the nuclear relaxation is or- 
ders of magnitude faster than previously thought. Our 
method is in general non-Markovian and is applicable to 
controls other than ultrafast lasers by appropriate choice 
of the Kraus operators. It can have wide application to 
other systems, such as gated QDs and NV centers in dia- 
mond 16 . An interesting application of the theory would 
be to use it for the design of the final nuclear state. 

This work was supported by LPS/NSA (EB) and in 
part by ONR and LPS/NSA (SEE). 
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